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^ Abstract 

We have performed single-atom magnetization curve (SAMC) measurements and inelastic scan- 

s 

ning tunneling spectroscopy (ISTS) on individual Fe atoms on a Cu(lll) surface. The SAMCs 
show a broad distribution of magnetic moments with 3.5 fiB being the mean value. ISTS reveals 
a magnetization excitation with a lifetime of 200 fsec which decreases by a factor of two upon ap- 
plication of a magnetic field of 12 T. The experimental observations are quantitatively explained 
by the decay of the magnetization excitation into Stoner modes of the itinerant electron system as 
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Magnetic atoms adsorbed on nonmagnetic substrates have been a topic of active study, 
both to provide insight into fundamental aspects of magnetism and as possible elements 
for future information technology and spin-based computation schemes. Depending on the 
type of substrate, ranging from thin insulating layers [IH3] over semiconductors p[|, to met- 
als [5-8J, one expects increasing hybridization of the atom with the substrate states. The 
itinerant nature of the substrate electrons play an increasingly pronounced role in the static 
and dynamic properties of the magnetic atom as one progresses through this sequence. 
Spin-polarized (SP-STS) jB] and inelastic (ISTS) [1] scanning tunneling spectroscopy or a 
combination of both J3J H] provide the only current means to probe magnetic properties and 
spin dynamics of isolated magnetic adsorbates on the atomic-scale. 

For magnetic atoms weakly hybridized with a non-magnetic surface, a description of the 
atomic moment by a half-integer spin governed by anisotropy terms in a spin Hamiltonian 
("isolated spin model") is sufficient [IHHIH]. Within this approximation, the role of the un- 
derlying host conduction electrons is neglected and the effect of the substrate is encompassed 
in the magnetic anisotropy. However, it has been known for decades that the magnetic mo- 
ments of 3d impurities in a nonmagnetic 3d, 4d or 5d metal are influenced qualitatively by 
the itinerant conduction electrons of the host material. Accordingly, as theoretically pre- 
dicted a long time ago, the magnetization dynamics of such systems is damped by decay 
into electron-hole pairs, namely Stoner excitations of the itinerant electron gas [10]. One 
consequence of this decay is a substantial g-shift and energy-dependent linewidth of the 
magnetization excitation if detected via a local method such as STS [TO] [EE]. Therefore 
the application of the isolated spin model to magnetic atoms on metals, while descriptive, 
inadequately describes these types of effects. This view is reinforced by the large excitation 
linewidth observed in previous ISTS which cannot be described by an approach which treats 
the magnetic atom as an isolated entity governed solely by magnetic anisotropy [7J IT2] . 
While the effects of Stoner excitations have been observed experimentally by spin-polarized 
electron-energy- loss spectroscopy of ferromagnetic surfaces [33] Q3], an experimental verifi- 
cation of their importance for the magnetization dynamics of individual impurities is thus 
far lacking. 

In this letter, we reveal the itinerant nature of individual Fe atoms absorbed on a Cu(lll) 
surface utilizing a combination of single-atom magnetization curve (SAMC) measurements 
[B] as probed by SPSTS and tunneling-electron driven excitations (ISTS). SAMCs reveal 
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the magnetic moment of individual Fe atoms which is ~ 3.5/ie. Complementary to this 
technique, ISTS reveals an anisotropy gap as well as a large linewidth of the magnetization 
excitation which increases linearly with magnetic field. In addition, we present the first 
theoretical studies which incorporate spin-orbit coupling in the response of an isolated local 
moment hybridized strongly with the host conduction electrons. The calculations provide a 
quantitative account of the anisotropy gap, and the linear variation with magnetic field of 
both linewidth and peak energy in the magnetization excitation spectrum. This confirms 
that the excitation lifetime is limited by decay into Stoner modes with a density of states 
linearly increasing in energy as predicted earlier (TUJ [TTJ . 

All experiments were performed on a home-built UHV STM with a base temperature 
of T = 0.3 K capable of applying a magnetic field B perpendicular to the surface [T5j . 
Non-spin sensitive W tips were used for ISTS and Cr-coated W tips were used for SAMC 
measurements [6j [16]. The Cu(lll) surface was cleaned by repeated Ar + sputtering and 
annealing cycles and Co islands were deposited at room temperature with a nominal coverage 
of 0.5 ML [TTJ . Such islands have single out-of-plane magnetized domains and were utilized 
to confirm the out-of-plane sensitivity of each Cr-tip [18]. Subsequently, Fe was deposited 
onto the cold surface at T < 6 K [19] . This results in a statistical distribution of single Fe 
atoms with a nominal coverage of a few thousandths of a monolayer (Fig. [Ija)). The atoms 
were observed to occupy only one adsorption site [20]. STM topographs were recorded in 
constant-current mode at a stabilization current / sta b with a bias voltage V^ ta b applied to the 
sample. SAMCs are recorded as described in Ref. [6] and in the supplementary material [21J. 
For ISTS, dI/dV(V) is recorded in open feedback mode via a lock-in technique with a small 
modulation voltages V mo d added to the bias (/ = 4.1kHz). 

In order to determine the magnetic moment of the individual Fe atoms, we measured 
the SAMCs of several single atoms (Fig. [jja)) using an out-of plane magnetized tip [6]. 
dl/dV maps acquired at various B of the same area reveal magnetic contrast on top of the 
Fe atoms [21]. From such maps, the spatially- resolved magnetic asymmetry is calculated 
by (dI/dV(B t ) - dl /dV{B x ))/{dI '/dV{Bf) + dI/dV(B±)) as given in Ref. [22]. Fig. [gb) 
illustrates a magnetic asymmetry map evaluated at opposite saturation fields (B^ = —0.4 
T, B± = +0.4 T). Each SAMC (e.g. Fig. [jjc)) is acquired by extracting dl/dV, for a given 
atom, during a B sweep (each B value corresponds to data extracted from one dl/dV map). 
Similar curves have been measured for about 60 different atoms at different locations with 
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different spin-polarized microtips where each atom probed had a minimum distance of 2 nm 
from any other magnetic structure. All reveal paramagnetic behavior with a saturation at 
B sat « 0.2 T. 

Because of the hybridization of the atom with the substrate, a quasiclassical continuum 
model j5j [6] is appropriate to describe the measured SAMCs. Fits utilizing various models 
produce a similar distribution of moments as a result of the rather significant out-of-plane 
anisotropy [23]. This model properly reproduces the measured curve if we assume an effective 
magnetic moment m ~ 3.5 /xb with a magnetic anisotropy of ~ 1 meV (Fig. [ijc)). The 
histogram of fitted m for all measured atoms is shown in Fig. [l](d). The average value is in 
good agreement with the DFT calculated total magnetic moment which considers both the 
spin and orbital moments of the Fe and the neighboring Cu atoms [231 124] . However, there is 
a broad distribution of the effective m. We speculate on two substrate-mediated effects: (i) 
there is a spatially varying mean field due to the inter-atomic long range RKKY interaction 
which mimics a different effective m for each atom [Ul [H] • (ii) The varying substrate density of 
states resulting from surface-state electron scattering might additionally change the value of 
m. While the detailed distribution of m changes when we use different minimum separation 
cutoffs, the maximum is always centered around 3.5 /xb independent of this choice. 

Complementary to SAMC measurements, ISTS can reveal information about the dy- 
namical magnetic properties of the Fe atoms. Fig. [2fa) shows ISTS curves taken on an 
isolated Fe atom (top) and on the Cu(lll) substrate (bottom) using the same tip atom as 
the magnetic field B is varied. While the substrate spectra are relatively flat and show no 
clear trend in a field, the atom spectra always show two steps symmetric with respect to 
the Fermi level Ep (V = OmV). These shift linearly towards higher voltage with increasing 
B and are accompanied by a simultaneous increase in linewidth. Numerical differentiation 
of the spectra illustrates both effects (Fig. |2](b)). Slight asymmetries in the lineshape for 
positive and negative voltage originate from small features in the tip or substrate density 
of states (Fig. ST a) bottom). Similar spectra using different tips were recorded on tens of 
different atoms with high energy resolution (« 150/xeV). All atoms show the same behavior, 
although with small quantitative differences in energy and linewidth. 

The peak energy E and full width half maximum (FWHM) linewidth AE of the step 
was extracted from Fig. |2](b) utilizing gaussian functions and plotted as a function of B in 
Fig. |2^c) and (d). As shown in the histogram in Fig. |2](c) (inset), the peak energy has an 



4 



average slope of 2.1 /ib (with a min/max of 2.1 ±1 /^b)- This proves that the observed steps in 
dl/dV are due to magnetic excitations of each Fe atom induced by tunneling electrons [U 12] 
resulting in a mean g-factor of 2.1 for this system. The height of the steps reveal that 
only about ~ 5% of all tunneling electrons induce such excitations |9J. The zero field 
excitation energy in the ISTS curves indicates that each Fe atom is subject to a substrate- 
induced magnetic anisotropy [U |2] of about 1 meV, which is four times smaller than the 
DFT predicted uniaxial out-of-plane anisotropy [23]. This measured anisotropy value was in 
turn utilized in the SAMC fitting model above. In passing we remark that, since the energy 
of the tunneling electrons used for the SAMCs is larger than the excitation threshold, 5% of 
the tunneling electrons excite the magnetization during recording the SAMC, leading to a 
time of 3 nsec between consecutive excitation events (7 s tab = 0.6 nA). This time is six orders 
of magnitude shorter than the time resolution of the experiment and more than four orders 
of magnitude longer than the lifetime of the excitation, justifying the assumption of thermal 
equilibrium within the model used for the SAMC curve fitting. 

Most notably, the excitation steps are very broad, with a width of AE = 1.5 meV that 
increases linearly with B as shown in Fig. |2^d). Since artificial experimental broadening 
due to Vmod and a nonzero temperature [25J is negligible in our experiment, the lifetime 
of the excitation is calculated to be remarkably short r = h/(2AE) = 200fsec (B = 0T). 
Comparably short lifetimes of magnetic excitations of single atoms have been reported for Co 
and Fe adsorbates on Pt(lll) and ascribed to the interaction between the atom states and 
the substrate electrons [TJ, [T5] . In these papers, the magnetic anisotropy and the relaxation 
mechanism were both described within an isolated spin model, neglecting itinerant effects, 
which is misleading in the case of strong hybridization. In contrast, both experiments 
reported here demonstrate the itinerant nature of the Fe/Cu(lll) system, as highlighted by 
the energy- dependent linewidth of the magnetic excitation and by the non-half integer total 
angular momentum J = m/{gyL-Q) = 1.75. Thus, a description of the Fe atom within an 
isolated spin model is not suitable. 

We turn next to a description of our theoretical studies, which incorporate the itinerant 
effects evident in the data. In our description of the magnetization excitations, we address 
the frequency dependent local transverse spin susceptibility The imaginary part of 

this function provides us with the density of states of magnetization excitations. Thus, 
Im[x(f2)] can be directly compared with the data displayed in Fig. |2](b). Two approaches 
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have been developed to compute \ for nanostructures: an empirical tight-binding approach 
(ETB) [TTj and a method [2"B"1 |2"7] based on the time-dependent density functional theory 
within the Korringa-Kohn-Rostoker (KKR) scheme [2E]- While the former requires param- 
eters, namely the various hopping integrals, the latter does not incorporate the challenging 
task of including the spin-orbit coupling which is essential to properly reproduce the ex- 
perimental observations. Thus, we report here our calculations obtained within the ETB 
approach taking full account of the spin-orbit interaction, which is a major advance from 
the theory presented in Ref. [UJ. Thus, the influence of magnetic anisotropy on the exci- 
tation spectrum is incorporated fully in our analysis, and the spin-orbit contribution to the 
g— shift is present as well. We have extended the formal development described in Ref. [29] 
to the description of the local susceptibility of a moment-bearing atom. We also used the 
KKR-based scheme described in Ref. [26J by mimicking the spin-orbit interaction with an 
external static magnetic field. Results for generated with such a method are in good 
accord with those obtained with the ETB scheme we present here. 

Before discussing the theoretical magnetization excitations, we note that the two methods: 
KKR and ETB, with parameters fitted to reproduce the electronic structure obtained with 
KKR, lead to ground state moments of 3.2 /zb and 3.6 /xb, respectively. 

In Fig. [3j we show the results of our theoretical studies, obtained with the ETB scheme, 
of the magnetization excitation spectrum of the Fe atom. In zero field, we see that the spin- 
orbit anisotropy renders the transition energy finite corresponding to a magnetic anisotropy 
of 1.1 me V which is very close to the experimental value (Fig. [2^c)). The spin-orbit coupling 
parameter used here is 34meV. This is 35% lower than utilized in Ref. [29], but compatible 
with the range of coupling constants found in the literature. Furthermore, use of a somewhat 
larger value of the spin-orbit coupling parameter, as discussed in Ref. [29], leads to an 
anisotropy of 2.5 meV which is too large as compared to the experimental result. 

The B dependence of the peak energy in the excitation spectrum (Fig. |3j[b)) agrees 
nicely with the measured data. The theory gives a g-factor of 1.8, slightly smaller than the 
experimental mean. The absolute value of the linewidth and its increase with the magnetic 
field are qualitatively reproduced in the theory (Fig.|3](c)), although with a smaller slope than 
in the measured data. A modest adjustment of the Fe/Cu hopping integrals would resolve 
this difference. The very large linewidths have their origin in the decay of the excitation 
into electron- hole pairs [TU]. Consequently, the theory provides a very good account of the 
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systematic features in the data. 

In summary, with a combination of SAMC and ISTS measurements, we reveal the itiner- 
ant nature of the magnetization excitation of an Fe atom on the Cu surface. The signature of 
itinerant effects is manifested by a non-half integer magnetic moment and a large, magnetic 
field dependent, linewidth of its excitation. Calculations of the magnetic field dependent 
local dynamic susceptibility reveal that the excitation lifetime is dictated by a decay into 
Stoner modes of the itinerant electron system whose density of states is linearly increasing 
with energy and thus with magnetic field. While Ref. [TU] did not incorporate spin-orbit 
coupling and addressed moments in the bulk of a paramagnetic metal, the data reported 
here confirms the itinerant picture set forth in this paper. This picture will apply for other 
magnetic atoms on 3d, 4d or 5d metal surfaces. [HI [261 [27] It has been predicted that the 
coupling to Stoner modes is less efficient for (i) weak overlap of the d-minority states of the 
atom with E-p and (ii) weak overlap of the d-majority states with the d-bands of the sub- 
strate [HI |2S] ■ Such considerations consequently can serve to predict material combinations 
with long lifetimes of the magnetization excitation which is important for future spintronic 
applications. 
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FIG. 1: (color online), (a) Topograph of a distribution of Fe atoms on Cu(lll). The colorbar 
spans 90 pm. (b) Spatially-resolved magnetic asymmetry map between saturation field values 
B = ±0.4 T. (I stab = 0.6 nA, V stab = -lOmV, V mod = 5mV (rms)). Color bar from a.u. to 0.025 
a.u. (c) Dots: SAMC extracted from the -B-dependent dl/dV maps by averaging over the signal 
from an individual atom. Solid line: fit to the data using the continuum model, (d) Histogram 
of effective magnetic moments extracted from the fit to the SAMCs of Fe atoms with minimum 
distance 2 nm to other atoms or Co islands. 
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FIG. 2: (color online), (a) ISTS at different I? as indicated (7 s t a b = InA, V sta b = — lOmV, 
Vmod = 0.1 mV (rms)). Top spectra were taken on an Fe atom. Bottom spectra were taken with 
the same tip on Cu (vertically shifted by — 14 nS). (b) Numerical differentiation of spectra taken 
on another Fe atom using a different tip (/stab = 2nA, Vs ta b = — lOmV, V mo( i = 0.05mV (rms)). 
(c) Energy and (d) averaged FWHM (high resolution data; error corresponds to the maximum 
and minimum value) of the excitation as a function of B extracted from gaussian fits to spectra 

(b) . The inset in (c) shows a histogram of the measured (7-factors determined from the slope of 

(c) divided by /^b- The values are extracted from spectra measured on 78 (c) and 16 (d) different 
atoms. 
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FIG. 3: (color online) (a) Calculated density of magnetization excitations for an Fe atom on 
Cu(lll) for various magnetic fields using the ETB approach. At zero field, spin-orbit induced 
anisotropy renders the excitation energy finite, (b) and (c) show magnetic field dependency of the 
peak energy and of the linewidth AE in the excitation spectrum (a). 
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